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A  plane  wave  emergent  from  a  scattering  layer  is  modeled  as  a  boundary 
value  problem.  The  wave  is  specified  at  z  by  a  unit  amplitude  and  a  phase 
prescribed  by  a  homogeneous  isotropic  random  function  with  a  power  law  spec¬ 
trum  and  Gaussian  statistics.  It  is  shown  that  the  statistics  of  the  wave  pro¬ 
pagating  in  free  space  away  from  the  boundary  cannot  be  Gaussian,  except 
perhaps  in  the  far  field  region  where  the  scattering  is  saturated.  It  is  also  shown 
that  the  amplitude  and  intensity  spectra  of  this  wave  each  satisfy  specific  scal- 
ing  relations.  Such  a  wave  has  been  called  a  diffractal  by  Berry  and  herein  we" 
further  develop^ome  of  his  ideas  on  diffractals,  in  particular  extending  his  dis¬ 
cussion  from  two  to  three  spatial  dimensions.  A  study  of  the  scintillation  index 
indicates  that  a  diffractal  is  not  an  unreasonable  model  of  a  radio  wave  passing 
through  the  ionosphere,  as  it  does  in  satellite  communication. 


1.  INTRODUCTION 


In  this  paper  we  investigate  the  statistical  properties  of  a  freely 
propagating  plane  wave  whose  phase  is  specified  on  a  transverse  plane  located  at 
z  =  0.  The  wave  is  propagating  in  the  z-direction  with  constant  frequency  u  and 
wavelength  \(=2n/k)  with  amplitude  if/(r,z,t)  =  0-<ulu(r.z)  such  that 

u(r,z)  =  at  z  =  Q+  (1.1) 

where  <p  (r)  is  a  function  of  only  the  transverse  coordinates  r  =  (x,y).  In  a 
great  many  geophysical  situations  the  phase  function  is  assumed  to  be  random 
so  as  to  describe  such  phenomena  as  the  twinkling  of  stars  due  to  atmospheric 
turbulence, 1  radio  wave  scintillation  in  satellite  communication  due  to  electron 
density  fluctuations  in  the  ionosphere,  interstellar  scintillation,  interplanetary 
scintillation  of  quasar  radio  sources  caused  by  electron  density  irregularities  in 

A 

the  solar  wind,  etc. 

The  observed  scintillation  or  fluctuation  in  the  amplitude  and  phase  of  the 
received  wave  results  from  interference  of  phase  points  along  the  wavefront  in 
the  z  =  0  plane  as  the  wave  propagates  away  from  the  boundary.  The  interfer¬ 
ence  pattern  is  determined  at  least  in  part  by  the  statistical  properties  of  <p{r). 
Booker  et  al.®  were  the  first  to  use  an  expression  of  the  form  (1.1)  to  model  the 
effect  of  a  localized  scattering  region  on  the  subsequent  propagation  of  a  scalar 
wave  away  from  the  scattering  region.  The  basic  assumption  they  employed  is 
that  the  individual  scatterings  that  a  wave  undergoes  as  it  propagates  through  a 
region  of  fluctuating  refractive  index  are  not  important;  only  the  net  change  in 
the  observables  are  significant.  The  physical  observables  far  from  the  scatter¬ 
ing  region  are  the  change  in  amplitude  and  phase  in  the  received  wave  relative 
to  a  reference  wave  that  has  not  been  so  scattered.  Thus  if  one  is  interested  in 
only  the  net  effect  of  the  scattering  region  on  the  received  signal,  the  entire 
interval  can  be  represented  by  a  thin  phase  screen  having  the  desired  bulk 


properties,  i.e.  the  phase  screen  gives  rise  to  the  same  overall  phase  shift  and 
amplitude  variation  as  the  physical  scattering  interval. 

Since  its  introduction  ( 1 .  l)  has  been  generalized  to  represent  an  extended 
scattering  region  by  building  up  a  sequence  of  thin  phase  screens.  °  Near  the 
screen  the  initial  field  fluctuations  affect  only  the  phase,  but  as  the  wave  pro¬ 
pagates  in  free  space  away  from  the  screen  an  interference  pattern  develops. 

The  distance  to  this  interference  pattern  is  dependent  on  the  density  gradients 
in  the  phase  screen.  The  larger  the  density  gradients  the  greater  is  the  angle 
through  which  the  transmitted  wave  is  scattered  and  the  nearer  to  the  screen  is 
the  subsequent  interference  pattern.  It  is  the  small  scale  spatial  structure  in 

Q 

the  scattering  medium  that  has  large  spatial  gradients  .  For  a  scattering  inter¬ 
val  whose  fluctuations  are  manifestations  of  some  turbulent  motion  there  is  a 
broad  spectrum  of  scales  associated  with  the  irregularities  in  the  medium.  For 
example,  the  spectrum  of  irregularities  of  the  electron  density  in  the  ionosphere 
has  a  range  of  from  a  few  meters  to  several  hundreds  of  kilometers,  i.e.,  five 
decades  or  more  in  spatial  scale.  Thus  radio  waves  with  frequencies  on  the 
order  of  10  MHz  and  lower  encounter  irregularities  on  the  same  scale  as  their 
wavelength  and  traditional  propagation  models,  such  as  geometric  optics,  can¬ 
not  describe  the  transmission  of  the  wave  through  the  phase  screen,  i.e.  multi¬ 
ple  scattering  effects  become  important.1®’11  It  is  remarkable  therefore  that 

weak  scattering  theory  has  been  so  successful  in  explaining  the  observed  scintil- 

1 2  13 

lation  of  radio  waves  in  satellite  communication  as  well  as  other  data.  ’ 

Model  calculations  based  on  phase  screen  descriptions  make  two  fundamen¬ 
tal  assumptions  about  the  fluctuations  in  the  medium.  The  first  is  that  the 
statistics  of  the  irregularities  are  described  by  a  zero-centered  Gaussian  distri¬ 
bution.  The  second  is  that  the  spectrum  of  scales  of  the  irregularities  has  a 
power-law  form.  These  two  assumptions  however  do  not  provide  a  complete 
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characterization  of  the  fluctuations.  One  must  also  introduce  an  outer  scale 
cut-off  in  order  to  avoid  a  long  wavelength  divergence  of  the  integrated  spec¬ 
trum.  Most  ionospheric  data  can  then  be  adequately  described  using  weak 

scatter  theory  with  a  power-law  spectrum  of  irregularities  with  power  exponents 

12  13 

between  3.8  and  4.8  and  an  outer  scale  cutoff  of  approximately  10  km. 

For  a  thin  screen  the  phase  in  (1.1)  is  usually  given  by  the  straightline  ray 

1 3 

path  integral  across  the  scattering  interval  (weak-scatter  theory.) 

a 

p(r)  =  f  N(t,z  )  dz  (1.2) 

-A£ 

where  N  (r ,z)  is  the  deviation  in  the  refractive  index  from  its  ambient  value, 

AZ,  is  the  thickness  of  the  scattering  interval  and  the  plane  wave  is  normally 
incident  on  the  ribbon  of  scatterers.  It  is  clear  from  (1.2)  that  the  statistics  of 
the  phase  fluctuations  are  the  same  as  those  in  the  index  of  refraction.  Thus 
<p(r)  is  a  Gaussian  process  only  if  the  irregularities  in  the  medium  are  also  Gaus¬ 
sian.  If  however  the  mechanism  leading  to  the  irregularities  in  the  medium  is 
turbulence,  then  as  is  well  known  the  statistics  of  the  fluctuations  cannot  be 
Gaussian.14  Therefore  either  <p{r)  is  not  a  Gaussian  process,  or  the  source  of 
the  power-law  behavior  of  the  spectrum  of  N-fluctuations  is  not  turbulence,  or 
both.  In  the  example  of  the  ionosphere  there  is  apparently  no  experimental  evi¬ 
dence  to  support  the  assumption  that  the  fluctuations  in  plasma  density  are 

normally  distributed,  although  in  situ  measurements  indicate  the  spectrum  is 

12 

a  power-law  over  a  restricted  range  of  scales. 

It  is  not  our  intention  here  to  provide  a  detailed  list  of  the  possible  instabil¬ 
ity  mechanisms  that  may  contribute  to  the  power-law  spectrum  of  density 
fluctuations  in  the  scattering  interval.  Rather  we  suggest  that  the  assumption  of 
Gaussian  phase  fluctuations  has  been  made  for  mathematical  convemence  and 
has  not  been  physically  justified,  at  least  in  the  case  of  weak  scattering  ^cf 
(1.2)].  In  spite  of  this  we  retain  (l.l)  to  represent  the  transmission  of  a  wave 
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through  a  region  of  tenuous  scatterers,  modeled  by  a  power-law  phase  screen 
and  assume  the  statistics  of  <p(r)  are  Gaussian.  We  show  that  the  power-law 
behavior  of  the  spectrum  implies  that  the  fluctuations  in  the  medium  satisfy  a 
scaling  relation.  We  assume  that  if  y  is  a  real  parameter  then  the  phase  function 
<p(yr)  can  be  obtained  from  p(r)  by  the  scaling  relation  <p{yr)  =  y4/l*<p{ r), 
where  d  is  the  Euclidean  dimension  (d  =  2  in  the  case  considered  here)  and  /z 
characterizes  the  scaling  behavior.  The  scaling  relation  is  a  geometric  con¬ 
straint  on  the  phase  surface. 

The  scaling  properties  of  waves  transmitted  through  fluctuating  media  with 

Gaussian  statistics  and  power-law  spectra  for  the  fluctuations  have  not  been 
2*5 

fully  explored.  Rino,  in  his  analysis  of  strong  scattering  in  a  power-law  phase 

screen,  investigated  the  results  for  the  second-order  statistics  of  intensity  using 

27 

a  model  developed  by  Gochelashvily  and  Shishov.  In  the  present  report  we 
also  obtain  results  for  second-order  statistics  of  intensity,  but  we  emphasize  the 
scaling  behavior  of  the  propagated  wave.  The  relation  between  the  present 
results  and  the  earlier  work  of  Rino  and  others  is  commented  on  in  the  appropri¬ 
ate  places  in  the  text. 

1*5 

Mandelbrot  has  recently  drawn  attention  to  the  physical  importance  of 
functions  having  scaling  properties  in  his  important  work  on  fractals.  We  apply 
some  of  his  observations  to  the  phase  surface  given  that  it  is  assumed  to  obey  a 
scaling  law.  One  consequence  of  the  scaling  (self-similarity)  property  of  the 
phase  function  <p{r)  is  that  irregularities  exist  on  all  scales  so  that  such  func¬ 
tions  do  not  possess  derivatives,  i.e.  one  cannot  construct  a  tangent  to  the  sur¬ 
face  p(r)  at  any  point  in  space.  Mandelbrot  not  only  advances  the  thesis  that 
most  physical  processes  are  discontinuous  and  described  by  such  non- 
differentiable  functions,  but  in  addition  those  processes  studied  by  Norbert 
Weiner,  i.e.  Gaussian  random  functions,  although  of  this  kind  are  in  fact  benign 


by  comparison.  Thus  although  Weiner  processes  are  presently  encountered  in 

practically  all  physical  models  of  phenomena  in  which  fluctuations  are  thought 

to  be  important,  Mandelbrot  contends  that  the  random  component  of  nature  is 

much  richer.  The  experimental  data  he  has  amassed  and  juxtaposed  to  support 

15 

this  contention  is  impressive. 

i  a 

Berry  has  applied  and  extended  some  of  the  ideas  of  Mandelbrot  in  the 
wave  propagation  context  by  introducing  the  notion  of  a  diflractal.  A  diffrac 
is  a  wave  that  has  encountered  a  fractal  object.  Such  an  object  has  a  self- 
similarity  property  such  that  it  possesses  structure  on  sill  scales.  An  incide  • 
wave  of  wavelength  X  is  fairly  insensitive  to  structure  on  scales  much  smaller 
than  X  and  for  structures  much  larger  than  X  the  techniques  of  geometric 
optics  can  be  employed  to  describe  the  wave-object  interactions.  However  this 
separation  of  effects  cannot  be  made  when  the  scattering  object  is  a  fractal 
since  then  the  self-similar  character  insures  that  there  will  be  structure  on  a 
range  of  scales  that  includes  X  and  thus  no  geometric  optics  limit  exists.  Berry 
models  a  diflractal  using  (1.1)  by  assuming  <p{x).  i.e.  <p  (r)  in  one  transverse 
dimension,  to  be  a  Gaussian  random  fractal  function  with  a  power-law  spectrum 
having  a  fractal  dimension  D  lying  between  1  and  2.  Physically  this  boundary 
condition  constitutes  a  thin  phase  screen  approximation  for  a  wave  reflected  by 
a  fractal  surface  or  refracted  by  a  slab  of  transparent  material  with  fractal 
refractive  index. 

In  Section  2  we  discuss  a  homogeneous,  isotropic  statistical  field  as  a  model 
for  the  fractal  phase  surface.  Some  of  the  physical  consequences  of  having  a 
power-law  spectrum,  which  one  does  for  a  fractal  function,  are  discussed  for  this 
case  including  the  relation  between  the  power-law  index  and  the  fractal  dimen¬ 
sionality  of  the  statistical  process.  In  Section  3  the  paraxial  approximation  is 
used  to  propagate  the  fractal  phase  surface  away  from  the  edge  of  the  phase 


screen.  The  statistical  properties  of  the  free  field  are  examined  by  calculating 
the  coherence  of  the  diffracted  wave  as  well  as  the  correlation  and  spectrum  of 
the  diffractal  intensity.  In  Section  4  the  scintillation  index  is  calculated  as  a 
function  of  distance  from  the  phase  screen  and  it  is  determined  that  a  diffracted 
cannot  have  Gaussian  statistics  except  perhaps  in  the  far  field  region  where  the 
fluctuations  are  saturated. 


I 


2.  HOMOGENEOUS  PHASE  SCREEN 

Herein  we  assume  that  the  phase  surface  ^(r)  is  a  locally  homogeneous 
random  field,  that  is  to  say  a  random  field  with  statistically  homogeneous  incre¬ 
ments  and  follow  the  discussion  of  Monin  and  Yaglom14  to  establish  its  general 
properties.  Spatial  homogeneity,  like  time  stationarity,  means  that  the  distribu¬ 
tion  of  the  increments  [?(rt)  -  ^(r2)  ]  depends  only  on  the  magnitude  of  separa¬ 
tion  between  the  two  phase  points,  i.e.  on  I  ri  -  r2  I .  Consequently  the  probabil¬ 
ity  density  is  invariant  to  any  rigid  translation  of  the  field.  It  is  easily  demon¬ 
strated  that  the  average  of  the  incremental  phase  [  p(r  +  R)  -  <p{ r)  ]  in  such  a 
case  is  a  linear  functional  of  R 


<  [  <p(r  +  R)  -p(r)  ]  >  =  c  R  (2.1) 

If  we  further  assume  that  the  field  is  locally  isotropic,  meaning  that  the  proba¬ 
bility  density  for  the  incremental  phase  is  invariant  under  rotations  and 
reflections  of  the  spatial  increments,  then  c  is  the  zero  vector.14  The  mean 
incremental  phase  is  therefore  zero. 

Inasmuch  as  we  are  also  assuming  that  the  incremental  phase  is  a  Gaus¬ 
sian  process  it  is  completely  determined  by  its  first  two  moments.  The  general 
second  moment  of  the  incremental  phase  is 


<  [  v  (r,  +  R,)  -  (r^  ]  [<?  (r2  +  R,)  -  ^(r2)  ]>  =  D  (r,  -  r2,  R„  Re) 

For  a  locally  homogeneous,  isotropic  phase  surface  (2.2)  can  be  written  as 


D{t\  -  r2,  R,.Re)  =  g- 


TX  -  r2  -  Ri  1 1  +  ri  -  na  +  Re  i  ] 


-  2?|l  Tj  -  r2  -  R!  +  1^  i  j  -  ri  -  r2l  j 


(2.2) 

14 

(2.3) 


where  D{R)  is  the  phase  structure  function 


D(R)  =  <[?(r+R)-?(r)]z> 


(2.4) 


The  phase  structure  function  is  as  important  for  the  incremental  process  as  the 
correlation  function,  p( R)  [  =  +  R)  y  (r)  ^  j,  is  for  the  phase  surface.  From 

the  definition  (2.4)  it  is  clear  that 

D{R)  =  2[p(0)-p(*)  ]  (2.5) 

so  that  if  we  know  the  correlation  function  for  the  process  we  know  the 

structure  function  as  well.  The  converse  of  this  statement  is  not  always  true 

however,  since  a  knowledge  of  Z?(R)  will  not  uniquely  determine  p  (R).  In  the 

situation  we  wish  to  consider  here  the  correlation  of  the  phase  surface  <p(r)  with 

p(r  +  K)  vanishes  as  R  -*  *  so  that 

Um  D(R)  =  DH  =  2p(0);  p(R)  =  [/?(«)-  D(R)\  .  (2.6) 

Thus,  for  locally  homogeneous  isotropic  random  processes  with  a  correlation 
function  p(R)  which  vanishes  at  infinity,  the  statistical  characteristics  of  p{R) 
and  Z?(R)  can  each  be  determined  from  the  other 

Let  us  now  consider  the  spectral  decomposition  of  the  phase  surface  p(r) 
with  homogeneous  isotropic  increments  and  the  corresponding  structure 
function  Z?(R): 

*(r)  =  /(e<k'r  -1  )<I4  (k)  +  •  (2-7) 

The  integral  is  over  the  entire  wavevector  space  since  d  <f>(k)  =  $(d k)  is  a  ran¬ 
dom  function  of  the  wavevector  k  defined  on  all  intervals  (  k,  k  +  cfk )  not  con¬ 
taining  the  origin  and  <p0  =  p(0)  is  a  constant  random  variable.  The  two 
dimensional  spectral  density  't'  (k)  of  the  phase  surface  p(r)  is  defined  by 

<  d$(k)  d  4*(k)  >  =  +  (k)  dzk  (2.8) 

where  't'(k)  is  a  non-negative  even  function  of  the  wavevector  k.  The  phase 
structure  function  (2.4)  can  be  expressed  in  terms  of  the  spectral  density 
using  (2.7)  and  (2.8)  as 


D(R)  =  2  f  [1  -cosk-R]  *(k)  dzk  .  (2.9) 

We  note  that  the  function  *(k)  can  rapidly  diverge  as  I  k!  -»  0  according  to 
(2.9).  However,  the  factor  (1  -  cos  k-R)  tends  to  zero  as  (k-R)2  so  that  as  long 


/  k2  'if  (k)  d2k  <  » 
lkf<*o 


(2- 10] 


the  phase  structure  function  remains  finite  for  finite  R. 

Herein  we  assume  that  the  phase  structure  function  is  given  by 

D{R)  =  ;  0  <:  fi  -s  2 


(2.n; 


where  CM  is  a  positive  constant  determining  the  strength  of  the  correlations 
and  fi  is  a  parameter  to  be  determined.  The  inverse  Fourier  transform  of  (2.9), 
with  the  phase  structure  function  (2.11),  yields 


*<*>  •  -5 i&r  ****•"" 


=  ^^-CMr(/x/2)r(l +/z/2)sin(/i7T/2)/ iKi^^  . 


(2.12; 


Thus  we  have  a  power-law  spectrum  for  the  incremental  phase  and  integrating 


(2.12)  over  an  unbounded  region,  i.e. 


<<e2>  =  f*{K)d2K 


(2.13; 


we  observe  that  the  mean  square  phase  surface,  <  ip2  >,  is  infinite.  The  phase 
surface  therefore  satisfies  the  conditions  for  being  a  Gaussian  fractal  function 
with  homogeneous  isotropic  increments  and  spectral  density  *(K). 

A  fractal  function  F(x)  is  a  continuous  but  non-differentiable  function  that 


satisfies  a  scaling  relation  of  the  form  F{yx)  =  yM/<1F(x)  where  y  is  a  real 
parameter,  d  is  the  Euclidean  dimensionality  of  the  space  and  yu.  is  a  real 
parameter.  Thus  in  the  case  d  =  1,  if  the  variations  in  the  function  are  known  in 


an  Interval  yxQ  a  z  &  x0,  they  can  be  determined  in  the  expanded  interval 
yz  z0  i  x  ^  y  z0,  as  well  as  in  the  contracted  interval  z0  i  x  &  xq/  y.  This  scaling 
behavior  is  known  as  self-similarity  and  is  one  of  the  defining  properties  of  frac¬ 
tal  processes.  A  process  described  by  a  function  the  graph  of  which  has  a 
HausdorflC-Besicovich  dimensionality  greater  than  unity  is  said  to  be  fractal,  i.e. 
a  fractal  process  is  a  process  described  by  a  function  that  exhibits  a  fractal 
dimensionality.  For  example,  the  fractal  dimensionality  D  of  a  Wiener  pro¬ 
cess  is  two  in  a  Euclidean  space  of  one,  two  and  three  dimensions.  In  general, 
fractal  functions  (processes)  with  D  <  2  and  D  <  d  are  not  space  filling.  In 
order  to  maintain  the  self-similarity  property  a  fractal  process  can  only  occupy 
space  in  clustered  or  localized  patches.  The  phase  function  $p(r)  is  such  a 
random  fractal  function  in  the  Euclidean  space  d  =  2.  Its  fractal  (Hausdorff- 
Besicovich)  dimensionality  is  shown  below  to  be  related  to  the  parameter  Jjl  by 
fi/2  -  2-D. 

The  strength  of  the  phase  structure  function,  i.e.  the  constant  CM  in  (2.12), 
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is  here  chosen  by  analogy  with  Berry.  We  note  that  the  vector  connecting  the 
end  points  of  the  phase  increment  [^(r  +  R)  -  $p(r)]  has  a  finite  slope,  even 
though  the  function  <p(r)  is  non-diflerentiable.  We  use  this  fact  to  define  the 
strength  of  the  fractal  as  the  distance  over  which  the  phase  structure  function 
has  a  slope  of  one  radian,  i.e. 

D{L)/LZ  =  1  .  (2.14) 

The  parameter  in  (2.11)  can  now  be  replaced  with  LZ_M  yielding  for  an  arbi¬ 
trary  separation  vector  R, 

D{R)  =  LZ*R»  .  (2.15) 

The  notion  of  a  fractal  (HausdorfT-Besicovich)  dimensionality  is  a  useful  con¬ 
cept  for  understanding  the  degree  of  irregularity  possessed  by  a  random  pro¬ 
cess.  We  use  a  simple  covering  argument  to  show  that  the  fractal  dimensionality 


D  of  the  phase  surface  is  related  to  the  exponent  m  in  (2.15)  and  (2.12)  by 
D  =  (  4  -  /x)/2.  The  argument  is  similar  in  spirit  to  one  presented  by  Mandel- 
brot,  wherein  he  uses  the  Lipschitz  condition  to  determine  the  fractal  dimen¬ 
sionality  of  a  given  process.  Instead  of  the  Lipschitz  condition  we  use  the  result 
for  the  phase  structure  function  (2.15).  We  note  that  the  number  of  spheres  of 
radius  R,  in  ($,r)  -space,  required  to  fill  the  volume  between  the  phase  surfaces 
at  r  and  r  +  R  is 

n  =  [D{R)]*RZ/  /?3~  R^z~l  as  R  -*  0  .  (2.16) 

The  total  number  of  spheres  required  to  fill  the  volume  between  r  =  0  and  r  =  1 
is  N  =  n/R  as  R  -*  0.  The  fractal  dimension  D  follows  from  the  definition 

N  ~  Rd  (2.17a) 

in  the  limit  R  -»  0  and  is 

D  -  log  N/  log  (1//?) 

=  (4-a*)/2  (2.17b) 

so  that  D  falls  in  the  interval  1  ^  D  «  2,  since  from  (2.11),  0  <>  fj.  ■&  2.  Berry 
obtained  the  same  result  in  the  one  dimensioned  case  using  a  capacity  argu¬ 
ment.  Thus  the  fractal  dimension  is  independent  of  the  Euclidean  dimensionality 
of  the  (  $,r)  -  space. 

As  pointed  out  above,  the  defining  property  of  a  fracted  function  is  its  self- 
similar  character.  This  is  also  true  of  a  random  fractal  function,  but  in  this 
latter  case  the  self-similarity  is  also  manifest  through  the  scaling  of  the 
corresponding  probability  density.  Let  P  (  $,r)  be  the  probability  density  that 
the  random  function  ^(r)  has  a  value  in  the  interval  (  $,  $  +  d$)  at  the  location 
r  .  Then  since  P  (  $,r)  is  Gaussian  and  ^(r)  has  the  phase  structure  function 
(2.15),  we  observe  for  a  positive  real  parameter  y  the  scaling  relation^0 


P($,r)  =  P{y)i/Z  $,  yr)  y^z 


(218) 


where  p.  is  the  parameter  in  (2.15).  Equation  (2.18)  indicates  that  if  the  phase 
surface  <p{r)  is  scaled  in  both  the  x  and  y  directions  by  7  and  in  the  $  direc¬ 
tion  by  7 then  the  resulting  phase  surface  is  statistically  indistinguishable 
from  the  original.  This  scaling  property  implies  that  in  addition  to  not  having  a 
smallest  scale,  which  is  required  for  a  fractal  function,  the  phase  function  p(r) 
also  has  no  largest  scale  and  is  self-correlated  over  an  arbitrarily  long  distance  [ 
cf.  (2.11)]. 14,16 

In  the  case  n  =  1  the  probability  density  (2. 18)  scales  in  the  same  way  as 

Brownian  motion;  the  phase  structure  function  is  proportional  to  R  [cf.  (2.11)] 

so  that  we  refer  to  D  =  1.5  as  the  Brownian  fractal  [cf.  (2.17b)].  The  case  n  =  0 

is  the  extreme  fractal  where  the  phase  surface  ^  is  on  the  verge  of  filling  a 

finite  volume.  Finally  the  case  p.  =  2  is  the  marginal  fractal,  where  the  phase 

is  "almost"  a  smooth  function.  If  instead  of  specifying  the  phase  structure 

function  (2.1)  we  had  instead  originally  assumed  a  spectrum  ot  the  form  (2.12), 

then  we  could  not  inverse  Fourier  transform  this  spectrum  to  obtain  (2.12)  for  p. 

=  2  because  the  integral  would  diverge.  Therefore  we  would  be  unable  to  choose 

the  coefficient  CM  as  done  in  (2. 15).  If  however,  we  can  identify  Cz  as  the  root 

mean  square  "slope”  of  displacements  of  any  length,  this  would  replace  the 
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length  L.  These  latter  comments  follow  closely  those  of  Berry. 
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3.  THREE  DIMENSIONAL  DIFFRACTALS 

The  propagation  of  a  diflractal  in  free  space  can  be  determined  by  solving 
the  three  dimensional  wave  equation  subject  to  the  boundary  condition 
u(r.z)  =  exp[ifc?(r)]  at  z  =  0*.  For  a  plane  wave  of  frequency  «  and 
wavenumber  k  the  wave  equation  reduces  to  the  Helmholtz  equation 

u(r,z)  =  0  (3.1) 

with  the  boundary  condition 

ti(r,  z=0+)  =  e<<s*W  ,  (3.2) 

and  V2  is  the  two  dimensional  Laplacian.  The  Helmholtz  equation  is  an  elliptic 
partial  differential  equation  so  that  in  order  to  find  the  wave  field  at  a  given 
point  in  space  one  must  solve  the  equation  for  the  field  at  all  points  in  space. 

One  often  approximates  the  Helmholtz  equation  by  a  parabolic  equation  having 
the  form  of  the  Schrodinger  equation  in  two  space  dimensions  with  the  time 
replaced  by  the  propagation  coordinate: 1 

u(r,z)  =  0  .  (3.3) 

The  boundary  condition  (3.2)  is  the  "initial  value"  of  the  solution  to  (3.3).  In  the 

parabolic  approximation  the  normal  to  the  phase  front  of  the  wave  is  assumed  to 
remain  close  to  the  direction  of  incidence,  i.e.  the  z-axis.  For  this  reason  the 
approximation  is  also  referred  to  as  the  paraxial  approximation.  The  condition 
for  the  validity  of  (3.3)  is  therefore  small  angle  scattering  and  can  often  be 
expressed  in  terms  of  a  correlation  length.  The  validity  of  (3.3)  will  be  discussed 
subsequently  since  its  applicability  to  describing  the  propagation  of  diffractals  is 
not  obvious. 
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3. 1  Average  Diffractal 

The  solution  to  the  parabolic  equation  (3.3)  subject  to  the  boundary  condi¬ 
tion  (3.2)  is  given  by  the  elementary  diffraction  integral 

,  ^  k  |j.  g»|g 

u(r,z)  =  e4**  f  dV  e  2*  et**(i*)  (3.4) 

Equation  (3.4)  describes  the  diffraction  of  the  phase  front  away  from  the  phase 
screen  point  (  r'.O)  to  the  observation  point  (r.z).  The  average  diffractal 
detected  at  (r.z)  is  determined  by  an  average  over  an  ensemble  of  realizations 
of  the  phase  irregularities: 

<u (r.z)>  =  rr—e**/  dV  e*  z*Jr  *'1*  <  et* *(*')>  .  (3 .5) 

27 rz  J 

However,  the  Gaussian  nature  of  the  phase  fluctuations  allows  us  to  write 

<  ai* *(!-)>  =  e  - **<**>/ 2  s  0  (3.6) 

where  the  last  equality  follows  from  the  fractal  nature  of  cp(r),  i.e.,  the  mean 
square  phase  variation  of  a  fractal  is  infinite  [  cf.  (2.12)  and  (2.13)]. 

3.2  Coherence  of  diffractala 

The  coherence  between  the  diffractal  wavefield  (3.4)  at  the  points  (r  +  R,z) 
and  (r.z)  is  given  by  the  coherence  function 

< /(R,z)>  s  <u(r  +  Rz)u’(r,z)> 

2ttz  j  j 

<  *(*")]>  .  (3  7) 

Using  the  conditions  of  Gaussianity,  isotropy  and  homogeneity  we  evaluate  the 
average  in  (3.7)  using  (2  4)  to  obtain  after  an  elementary  transformation  of 


variables 
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* 


{ 


</(  Rz)> 


/  d*-q  f  dUe'k*D{i)/2 


(3.8) 


-  exp 


(R-f)  +  i~K  (R-e) 


The  integral  over  17  yields  (2 n-)2  <5  (R  -  f)  za/A:2  so  that  using  (2.15),  (3.8) 
reduces  to 


</(K)>  =  exp 


- k 2 


(3.9) 


which  is  independent  of  z.  The  coherence  function  (3.9)  is  consistent  with  the 

5 

result  first  obtained  by  Booker  et.  al.  for  a  thin  phase  screen  with  Gaussian 

statistics.  Most  recently  Berry  has  obtained  the  one  dimensional  form  of  (3.9) 
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for  a  diffractai. 

The  energy  at  each  scale  in  the  diffractai  is  determined  by  the  power  spec¬ 
trum  Pu{YCj  and  is  obtained  by  taking  the  Fourier  transform  of  the  coherence 
function  (3.9): 


PU(K)  =  yd*ue1E'*ex  p 


-i2/,2"**  luk/2 


(3.10) 


Note  that  (3.10)  has  the  form  of  a  radially  symmetric  Ldvy  stable  probability 

density  in  two  dimensions,  that  is  to  say  that  (3.9)  has  the  form  of  a  Ldvy 
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characteristic  function  with  parameter  p.  A  Ldvy  stable  process  is  one 

whose  distribution  is  a  member  of  the  class  of  infinitely  divisible  distributions 
and  therefore  it  satisfies  a  number  of  scaling  relations.  The  identification  of  the 
power  spectrum  of  the  diffractai  with  such  a  process  is  not  accidental;  rather  it 
emphasizes  the  lack  of  a  fundamental  scale  in  the  diffractai  spectrum.  The  scal¬ 
ing  properties  of  the  spectrum  are  made  evident  by  scaling  the  variables  in  the 
integrand  of  (3. 10)  by  [kz  L2'M/  2]1/M  and  integrating  over  angle  to  obtain 


PU(K)  =  2rr 


k 2  Lz~* 


f  y  dy  e  v*  J0 
a 


kzL2-» 


1/M 


'Kiy 


(3.11) 
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where  /0( )  is  the  zero  under  Bessel  function. 

In  the  limit  K  -*  0  the  Bessel  function  in  the  integrand  of  (3.11)  is  replaced 
by  unity  and  we  obtain,  in  terms  of  the  fractal  dimensionality  (2. 17b), 

2n  V  [(3-D)/ (Z-D)]  .  (3.12) 

Note  that  (3.12)  is  independent  of  spatial  wavenumber  and  that  the  spectral 
level  for  large  scale  features  scales  as  The  fractal  dimension  D  =  2  is 

somewhat  special  and  will  be  discussed  in  Section  4.  In  the  short  scale  asymp¬ 
totic  limit,  K  -»  «,  the  well  known  asymptotic  expansion  for  a  stable  Ldvy  pro- 
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cess  in  two  dimensions  applied  to  (3.11)  yields  ' 

]^mk2Pu  (K/fc)  =  2?-*°  (kL)W-"r2(3-D)  sin(2-Z7)n/ Ik/ (3.13) 

For  intermediate  values  of  the  wave  number  a  more  detailed  analysis  of  the 
integral  (3.11)  is  required.22'23 

The  asymptotic  results  for  the  diflractal  spectrum  (3. 12)  and  (3. 13)  can  be 
used  as  an  a  posteriori  justification  for  the  paraxial  approximation.  The  para¬ 
bolic  equation  is  an  adequate  description  of  the  propagation  away  from  the 
phase  screen  where  the  scattering  angle,  as  measured  by  !  Kl  /k,  is  small. 
Therefore  we  require  that  most  of  the  energy  in  the  scattered  wave  field  is  in  the 
forward  direction.  More  quantitatively,  if  I K  /k!  is  large,  we  require  that 

Pu(0)»P»  (  K/k): 

IK/*I-«1.  (3.14, 

Therefore  we  introduce  the  parameter  0  given  by 

2-m  D-\ 

0  =  (kL)  »  =  ( kL)2~D  (3.15) 

and  the  paraxial  approximation  is  a  good  one  if  k  9/  IkI«1,  or  more  strongly  if 
9«l.  Thus  for  a  given  incident  wavelength  X  for  L  sufficiently  small  the 


lim  kzPu(K/  k)  = 
K-0 
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dominant  diffracted  waves  are  paraxial.  The  scattering  parameter  0  increases 
with  increasing  wavenumber  (decreasing  wavelength)  of  the  incident  wave  indi¬ 
cating  the  small  scale  divergence  of  the  slope  of  the  phase  surface. 


3.3  Coherence  of  diffractal  intensity 

The  coherence  between  the  diffractal  intensities  I  u  (r  +  Rz )  1 2  and 
lii(r.z)l2  is  given  by 


</2(Rz)>  = 


2  n  z 


/  n  d*Ti 

J  >=i  1 


exp 


i  ~  [lr+R-r1l2-lr+R-r2lz  +  lr-r3l2-lr-r4lz] 


•  <  expi&  r,)  -*(**)  +  «»  (r3)  -  ?(r4) 


(3.16) 


The  average  in  (3.16)  requires  the  specification  of  the  mean  square  four  point 
difference  in  the  phase  surface 

<  [?(ri)  -?(r2)  +  <p{r£)  -?(r4)]Z>  =  Z?(lri  -  r;J)  +  £){\rl  -  r4l)  +  £(!r3  -  r4l) 


+  Z7(lr3  -  r2l)  -  Z7(lri  -  r3l)  -  D{\tz  -  r4l) 


(3.17) 


so  that  using  the  local  homogeneity  and  isotropy  assumption  on  the  phase  struc¬ 
ture  function  [cf.  (2.15)]  we  obtain 


</z(Rz)>  = 


k 

2nz 


f  it  dh-i 

J  j  =  i  1 


exp 


i  ^-|lr  +  R*  ri!z  -lr  +  R-r2!z  +  lr-r3!z  -lr  -  r4!2 1 


exp 


-  ^-kz  LZ'M  |ln  -  r2lM  +  lr,  -  r4!M  +  lr3  -  r4i“ 


+  ir3  - - i r i  -r3iM  -  ir2  - r4i' 


(3.18) 


The  linear  transformation  of  variables: 


*l\  =  |-(rt  +  r2  +  r3  +  r4).  rfz  =  |-(r,  -  r2  -  r3  +  r4). 

Va  =  +  r2  -  r3  -r4)  and  tj4  =  |"(rt  -  r2  +  r3  -  r4) 

enables  us  to  directly  integrate  (3.18)  over  iji  to  obtain  {  2n)2  6(rj4)  z2/ k2  so 
that  integrating  over  rj4  (3.18)  reduces  to 

</2(Rz)>  =  /  d2rjz  f  e  '**  **  ^ 

k  i 

•exp  j -i-ifczL2-^  [2T7^  +  2rj^-ln2+^3lM- (3.19) 

A  more  compact  form  of  (3. 19)  is  obtained  by  introducing  the  dimensionless 
parameter 

<■3  te9z  (3.20) 

and  using  the  dimensionless  variables 

u  3  k  Qrjz/t*.  v  =  k  0T73/  (3.21) 

we  obtain 

<  /2(Rz)  >  =  f  d*uf  d*v 

•  exp  iu- [v-Rfc  -fM/2/(u.v)  (3.22) 

where 

f  (u,  v)  =  lui^  +  W*-  J^u-vi^  .  (3.23) 

The  spectral  density  of  the  diffractal  intensity  in  terms  of  the  Fourier  transform 
of  (3.22)  can  be  written 

Pi  (K  0  =  f  d2R  eiI  B  <  [lit(r  +  Rz)  I2  -  <  iu(r  +  Rz)l2  >  ] 

•  f  I  u  (r,  z )  1 2  -  <  lu  (r,  z )  !z  >  ]  >  j 
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=  /  d*R  eirB  <  /2(R<-)  >  -  1 


(3.24) 


where  we  have  used  (3.9),  i.e.  <lul2>  =  1.  In  terms  of  the  scaled  wavevector 
q  =  YL/kQwe  can.  after  integrating  over  R  rewrite  (3.24)  as 

kzP,{ q*9,f)  =  -p-  /  d2!/  cos(v'  q)  [«-/(«**>  -  1  .  (3.25) 


The  spectrum  of  intensity  fluctuations  therefore  has  the  following  limiting  forms 
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obtained  by  Gocheiashvily  and  Shishov. 

In  the  low  frequency  limit  I  Kl  /  k  0  <  f  -1  ,  the  exponential  in  (3.25)  can  be 
expanded  and  the  first  order  correction  to  the  spectrum  given  by 

lim*2P/(qA:9,fl  a  ^-T2  (1  +  (jS  2)  lqlz  <«  \  q  <  j-  (3.26) 
q-a  nvr  C 

Thus  the  spectrum  goes  to  zero  quadratically  in  the  parameter  f .  and  as  the 
(2-/i)  power  of  the  wavenumber. 

In  the  high  frequency  limit  the  intensity  spectrum  can  be  evaluated  by  seg¬ 
menting  the  integral  (3.25)  into  a  part  v  <.  f  q  and  a  part  v>fq: 

kz  Pi  (t if)  =  ^rfdv  \  f  cos  (v  q)  e 
"o  l  o 

m  m  ' 

+  f  vdv  cos  (▼•  q)  e  'ffad  -  /  vdv  cos(vq)  r  (3.27) 

hr  o 

In  the  first  integral  we  expand  f  (  v.^q )  for  £  ■+  *>  and  obtain 

Um  / (v.Cq)  ~ 

in  the  second  integral  since  v  >  f  q ,  we  use 


Thus  (3.27)  yields 


lim  /(▼.<" g  )  ~  0  . 
C-o 


limfc2P/(q,t)  =  u(q/2l/M) 


v  *-*  r  *'•  •'  V*'-  Vs/  ‘  *•  *\*  \  *\  >  ' 

m*  _ **«.  »*_  ■*"  '  *“  .*  *  ■  •*  *  •  '  •  •  v  ■ 
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=  -\—Pu(  q/24^)  (3.28) 

22~d 

so  that  the  intensity  spectrum  saturates  far  from  the  phase  screen  resulting  in 
a  scaled  version  of  the  spectrum  of  u  [cf.  (3.14)]. 

As  a  final  comment  on  the  intensity  spectrum  we  observe  that  (3.25)  itself 
satisfies  a  scaling  relation.  Recalling  that  q  =  K/  k9  in  (3.25)  and  defining 

fczP/(K/fc.0,O  s  /  §-cos  (▼•  K/kQ)  fe  /(*.<•**«)  -  ij  (3.29) 

0 C  l 

we  obtain  the  scaling  law 

Pj(K/ k ,®,{)  =  P/(K//fc0,l,^)/02  .  (3.30) 

Thus  as  the  phase  surface  deformation  increases,  as  measured  by  the  parame¬ 
ter  0  ,  spectral  information  appears  closer  to  the  phase  screen  and  farther 
from  the  z-axis,  i.e.  at  larger  angles. 
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4.  SCINTILLATION  INDEX 

A  parameter  often  used  to  measure  the  strength  of  the  intensity 
fluctuations  is  the  scintillation  index  S4  defined  by  the  integrated  intensity 
spectrum: 

Si  =  -fPjiKOd^K  (4.1) 

Substituting  (3.24)  into  (4.1)  we  obtain  an  expression  for  the  scintillation  index 
in  terms  of  the  fourth  moment  of  u,  i.e. 

Si  =  <  /Z(R  =  Q,z)>  -  1  .  (4.2) 

Rumsey24  introduced  an  intensity  randomization  factor  U,  obtained  from  (4.1) 
by  replacing  f7(K,f)  with 

Pj  (Kfl  =  4Jfc2  LL  sin2  (K**/2k)  $„(K)  (4.3) 

where  A  L  is  the  thickness  of  the  phase  screen  and  (X)  is  the  spectrum  of 
irregularities  in  the  screen.  Introducing  a  turbulence  strength  parameter  T  by 
$  n  (K)  =  2  TtT/[k2&L  IkIm  f  2  ]  the  intensity  randomization  factor  can  be  writ¬ 
ten 

pZ-fi/Z  m 

U  '  (fcajr  T  sin2  {£qz) 

2  t*  1M 

=  T  —  T(1-m/2)  cos  (/xrr/4)  (4.4) 

in  agreement  with  Rumsey  since  (<$/  A:©)**  =  (z/fc)^2  In  his  analysis  of 

Ot N 

strong  scattering  from  a  power-law  phase  screen  Rino  determined  that 
Si  >  U  for  all  values  of  p.  except  in  a  small  region  around  p.-Z.  (This 
corresponds  to  his  parameter  value  v  =  1.5.)  For  the  fractal 
dimension  in  the  interval  1  <  D  ^  2,  the  v  parameter  is  confined  to  the  interval 
1.5,  i.e.  u  =  2.5-D.  At  the  end  of  the  v  interval  where  U  >  S2 ,  the 
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fractal  dimensionality  is  marginal,  i.e.  D  =  1.  The  marginal  fractal  separates 
diffractals  from  ordinary  random  waves.  The  physical  properties  of  these  two 
types  of  waves  are  quite  distinct  as  indicated  by  the  discontinuity  in  the  scintil- 

pc 

lation  index  obtained  by  Rino.  This  corresponds  to,u  =  2  in  (3. 19)  and  results 
in  </2(Rz)>  =  1,  so  that  the  saturation  level  </2(r,z)>  =  2  is  only  obtained  if 
the  £  -»  <*>  limit  is  taken  before  the  /i,  -»2  limit  [cf.  (4. 12)  below]. 

In  the  present  formulation  the  scintillation  index  (4.2)  can  be  written,  using 
(3.22)  with  R  =  0  and  the  symmetry  properties  of  the  integral,  as 

5f  =  ■— rg-  Re  f  dtu  fd?v  exp  iu  v  -  f*''2  /  (u,v)  -1  .  (4.5) 

In  the  region  close  to  the  edge  of  the  phase  screen  (  f  -*  0)  we  can  expand  the 
exponential  term  in  the  integrand,  of  (4.5)  to  obtain, 

Si  =  -^-Re  f  dhi  f  d*v  e4*  T 
4  (2jt)2  j  j 

•  1  [2u>*-lu  +  vl*]  +  0(?*)-l  .  (4.6) 

TTie  first  term  inside  the  curly  brackets  yields  (2rr)  2  <5  (u)  for  the  v  -  integra¬ 
tion.  The  u-integ ration  of  this  term  yields  unity  exactly  cancelling  the  -1  factor. 
The  term  involving  u**  yields  a  5  (u)  for  the  ^integration  and  <5(u)  vanishes 
for  the  u-integration.  The  integral  involving  >u  +  vlM  can  be  transformed  using 
x  =  )£{u  -  v)  and  y  =  )$(u  +  v)  to  obtain 

S l  =  ^C'Re  /  e  f  dzx  \x\»eix*  +  O(^) 

(2  it)4  j  j 

-  2?  «■  I  ^2  r  ( 1  +  2)  C0S(MTT/  4)  +  0{{»)  .  (4.7) 

From  the  power  spectral  density  (2.12)  we  observe  that  Rumsey's  turbulence 
intensity  factor  can  be  written^4 


o 
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t  =  r(/x/ 2)  r (1+^/2)  sin(^rTr-2) 


=  a»-‘M*ey  <“a> 

which  when  substituted  into  (4.4)  yields  one-half  of  (4.7).  Thus  we  observe  that 
Sf  =  2 U  in  the  near  field  region,  i.e.  the  region  close  to  the  phase  screen. 

If  instead  of  the  parameter  L  introduced  in  (2.14).  we  had  instead  selected 
CM  in  (2.12)  in  such  a  way  that  the  spectral  strength  is  given  by  Cp  as  defined 
by  (5)  in  Rino^9,  then 


and 


(4.9a) 


C  =  rq-M/2)  = 

M  n  +  fx/2) 


(4.9b) 


Substituting  (4.9a)  into  (4.4)  and  (4.9b)  into  (4.7)  again  yields  S?  =  2U  ,  in 
essential  agreement  with  Rino's  result  away  from  the  marginal  fractal  region 
(D  =  1  ). 


Finally  we  can  express  (4.7)  in  terms  of  the  fractal  dimension  as  follows: 


Si  =  26~2D  f~D  cos  T(3-Z7)  +  0(^2J})  .  (4.10) 

Equation  (4.10)  agrees  with  the  one-dimensional  result  of  Berry  which  for  the 
marginal  diffractal  D=l  +  s,e-»0,  yields 

5?  =  4t7£(-  +  0{e .  (4.11) 

The  scintillation  index  therefore  has  an  infinitesimally  slow  linear  growth  in 

£  (=kz)  for  D  -  1+s.1® 

Far  from  the  phase  screen  ((■-*»)  the  behavior  of  the  scintillation  index 
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is  quite  different.  Gochelashvily  and  Shishov  find  the  scintillation  index  to  be 
given  by 
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Q 


S$  =  1  +  C4(fj)  +  o(rz(z ■<*>) 

where 


(4.12) 


c4(/i)  s  2M r( i  +  m/2)  r ( 4//x - 1) 


1 

/MT 


r(  1  +  /</  2)  sin  /tTT/  2  + 


2  -  u 

4f(2-/i/2) 


•  [2r(/i/2)r(l-/i/2)]2^-z)/*‘  .  (4.13) 

Thus  the  scintillation  index  for  the  fractal  surface  approaches  unity  from  above 
as  f  [cf.  (4.12)3  so  that  the  intensity  scintillations  become  saturated  with 

propagation  distance  and/or  with  the  increased  strength  of  "turbulence" 

[cf.  (4.9)]. 

We  now  use  the  scintillation  index  to  establish  that  a  diffractal  caranot  have 
Gaussian  statistics  except  perhaps  in  the  region  far  from  the  phase  screen  To 
establish  this  we  make  the  contrary  assumption,  i.e.  that  a  diffractal  can  be 
described  by  a  complex  Gaussian  wave  field,  and  demonstrate  that  this  assump¬ 
tion  leads  to  a  value  of  Si  that  is  inconsistent  with  (4.10)  and  (4.12).  Thus  we 
write  the  diffractral  wave  amplitude  in  terms  of  its  real  and  imaginary  parts 

u(r,z)  =  ^(r.z)  +  uLi(r,z)  (4.14) 

where  Vy.  and  u*  are  real  homogeneous  Gaussian  random  variables.  Then  the 
mean  square  amplitude  is 

<itz>  =  <u?>  -  <U(Z>  +  Zi<urui>  (4  15) 

where  the  bracket  here  denotes  an  average  over  the  Gaussian  distribution.  But 
we  know  that  <  ua  >  =  0  because  of  the  fractal  nature  of  the  phase  [cf.  (3.6)],  so 
that  the  components  of  the  field  are  statistically  independent 

<urui>  =  0  (4.16a) 

and  their  mean  square  levels  are  equal 


w 


<u?>  =  <u?>  =  %  (4.16b) 

where  the  value  stems  from  the  normalization  (3.9).  Now  we  can  write  the 

fourth  moment  as 

sf  +  i  =  <U*> 

=  <  Ur  *>  +  <Ui*>  +  2<UrZ  Ui>  (4. 17) 

and  using  the  well  known  properties  of  Gaussian  functions14-  we  obtain  for  the 
scintillation  index 

S|  =  3<iir2>2  +  3<ui2>2  +  2<,uTz><iii2>  +  4<iiiur>2  -  1  .  (4.18) 

Inserting  the  second  moment  values  from  (4.16)  into  (4.18)  we  obtain 

Si  =  1  (4.19) 

for  a  homogeneous  isotropic  Gaussian  wave  field. 

However  the  scintillation  index  for  the  diffractal  given  by  (4.12)  is 


Si  -  1  + 


A 

2z(a-i> 


1 

z<D~l) 


(4.20) 


where  the  definition  of  the  seeding  parameter  (  has  been  used  [cf.  (3.20)]  and  A 
is  a  z  independent  positive  definite  constant.  Since  U  Ds2  the  scintillation 
index  (4.20)  approaches  that  for  a  Gaussian  random  wave  as  z  -»  <*>.  This  implies 
that  the  statistics  of  a  diffractal  are  not  Gaussian  except  perhaps  in  the  far 
field  region  where  the  fluctuations  are  saturated. 


5.  SUMMARY  AND  CONCLUSION 
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In  an  earlier  report  we  have  argued  how  the  existence  of  different  scales 
of  fluctuations  in  the  ionosphere  can  give  rise  to  the  observed  power-law  spec¬ 
trum  of  plasma  irregularities.  This  property  of  the  ionosphere  suggests  that  it 
can  be  reasonably  well  modeled  as  a  region  with  a  fractal  index  of  refraction. 
Subsequently  a  wave  transmitted  through  such  a  medium  would  become  a 
diffractal  as  proposed  by  Berry.  We  have  shown  herein  that  a  diffractal  results 
from  the  diffractive  free  space  propagation  away  from  a  boundary  of  a  homo¬ 
geneous.  isotropic  statistical  phase  surface  having  a  power-law  spectrum.  The 
index  /x  of  the  power-law  spectrum  at  the  boundary  is  related  to  the  fractal 
dimensionality  (D)  of  the  diffractal  by  D  =  (4-  /x)/ 2.  The  statistical  properties  of 
the  free  field  were  determined  by  calculating  the  coherence  of  the  diffractal  as 
well  as  the  correlation  function  and  spectrum  of  the  diffractal  intensity.  These 
low  order  moments  indicate  that  the  diffractal  can  be  characterized  by  scaling 
laws  in  the  spectrum  of  the  wave  amplitude  and  in  the  intensity  spectrum. 

Finally  the  scintillation  index  Sf ,  often  used  as  a  measure  of  the  strength 

of  the  intensity  fluctuations,  is  expressed  in  terms  of  the  distance  from  the 

phase  screen  boundary.  For  a  wave  field  with  Gaussian  statistics  Sf  =  1 

independent  of  z.  and  we  have  shown  that  a  diffractal  approaches  this  value  from 

above  as  an  inverse  power  of  the  distance  from  the  phase  screen  in  agreement 
2*5 

with  Rino.  This  implies  that  the  statistics  of  a  diffracted  are  not  Gaussian 
except  perhaps  in  the  far-field  region. 
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